A Shell Model Description of the Decay Out of the Super-Deformed Band of 36 Ar 
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Large scale shell model calculations in the valence space spanned by two major oscillator shells 
(sd and pf) describe simultaneously the super-deformed excited band of 36 Ar and its spherical 
ground state. We explain the appearance of this super-deformed band at low excitation energy 
as a consequence of the very large quadrupole correlation energy of the configurations with many 
particles and many holes (np-nh) relative to the normal filling of the spherical mean field orbits 
(Op-Oh). We study the mechanism of mixing between the different configurations, to understand 
why the super-deformed band survives and how it finally decays into the low-lying spherical states 
via the indirect mixing of the Op-Oh and 4p-4h configurations. 

PACS numbers: 21.10.Sf, 21.60.Cs, 23.20.Lv, 27.40.+Z, 29.30.-h 



Excited deformed bands in spherical nuclei provide a 
spectacular example of coexistence of very different struc- 
tures at the same energy scale, that is a rather peculiar 
aspect of the dynamics of the atomic nucleus. Several 
cases are known since long, for instance the four particle 
four holes and eight particles eight holes states in 16 O, 
starting at 6.05 MeV and 16.75 MeV of excitation energy 
0,0- However, it is only recently that similar bands, de- 
formed and even super-deformed, have been discovered 
in other medium-light nuclei such as 56 Ni 36 Ar p| 
and 40 Ca 5] and explored up to high spin. These ex- 
periments have been possible thanks to the advent of 
large arrays of 7 detectors, like Gammasphere or Eu- 
roball. One characteristic feature of these bands is that 
they belong to rather well defined spherical shell model 
configurations; for instance, the deformed excited band 
in 56 Ni can be associated with the configuration (lf 7 / 2 ) 12 
(2P3/2, lf 5 / 2 , 2P1/2) 4 while the (super)deformed bands 
in 36 Ar and 40 Ca are based in structures (sd) 16 (pf) 4 
and (sd) 16 (pf) 8 respectively. The location of the np- 
nh states in 40 Ca was studied in the Hartree-Fock ap- 
proximation with blocked particles and Skyrme forces in 
ref. 0- While many approaches are available for the 
microscopic description of these bands (Cranked Nilsson- 
Strutinsky Q , Hartree-Fock plus BCS with configuration 
mixing p|, Angular Momentum Projected Generator Co- 
ordinate Method @|, Projected Shell Model 0], Cluster 
models 0] e t c the interacting shell model is, when af- 
fordable, the prime choice. The mean field description 
of N=Z nuclei, has problems related to the proper treat- 
ment of the proton-neutron pairing in its isovector and 
isoscalar channels. On the shell model side, the prob- 
lems come from the size of the valence spaces needed to 
accommodate the np-nh configurations. 

In this letter we focus in the 36 Ar case, where a ro- 
tational band has been experimentally found starting 
4.3 MeV above the spherical ground state. It is gener- 
ated by an intrinsic super-deformed state i. e. bearing an 
axis ratio 2:1. In previous works [4I lllj. we have shown 



that the experimental super-deformed (SD) band can be 
associated to the promotion of four particles across the 
Fermi level, from the sd to the pf shell. Indeed, this 
is a very crude description, because the physical states 
contain components belonging to other np-nh configura- 
tions. This mixing allows for transitions connecting the 
superdeformed band and the low-lying sd states. How- 
ever, the mixing must be gentle enough so as not to jeop- 
ardize the very existence of the band. This is the theo- 
retical challenge that we affront in this letter: to explain 
why the super-deformed band appears at such low exci- 
tation energy and how can, at the same time, mix with 
other configurations and preserve its identity. 

Let's recall the results of the calculations at fixed 4p-4h 
configuration. In the description of the SD band of 36 Ar, 
the natural valence space comprises the sd and the pf 
shells. However, the inclusion of the ld 5 / 2 orbit produces 
a huge increase in the size of the basis and massive center 
of mass effects, forcing us to exclude it from the valence 
space. The quadrupole coherence of the solutions will 
be reduced by this truncation. We have quantified the 
effect and it is moderate as we shall show below. There- 
fore, our valence space consists of the orbits 2si/2, ld3/2, 
H7/2, 2p3/2, lfg/2j ancl 2pi/2- The basis dimensions are 
O(10 7 ). The starting effective interaction is the same 
used in ref. Q dubbed SDPF.SM and described in detail 
in 0]. Blocking the orbit ld 5 / 2 minimizes the center of 
mass contamination of the solutions. Nevertheless, to be 
safer, we have added to the interaction the center of mass 
hamiltonian with a small coefficient Xr m=0-5 to reduce 
even more unwanted mixings (see |l3j for an updated 
discussion of the center of mass issues) . The differences 
between the results with and without the addition of the 
center of mass hamiltonian are minor. 

In Figure n the calculated transition energies are com- 
pared with the experimental results in a backbending 
plot. The agreement is really remarkably good, except at 
J=12 where the data show a clear backbending while the 
calculation produces a much smoother upbending pat- 
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FIG. 1: The superdeformed band in 36 Ar; E 7 's, exp. vs. 4p- 
4h calculation 
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FIG. 2: The superdeformed band in 36 Ar; B(E2)'s, exp. vs. 
4p-4h calculation. The J=16 point is only a lower bound 



tern. In the experimental data there is a close-by sec- 
ond 10 + state, therefore, the discrepancy may be due 
to the lack of mixing in the calculation. If we move 
now to Fig. |3 where we gather the experimental fll| 
and calculated B(E2)'s (we use standard effective charges 
Sq^=5q u =0.5 and b=1.94 fm), we find again an aston- 
ishing accord The value of the intrinsic quadrupole 
moment corresponds roughly to a deformation (3=0.5, so 
that it is justified to speak of a super-deformed band. 
The origin of such a strong coherence can be attributed 
to the dominance of Elliott's SU(3) structures in the va- 
lence space. The almost perfect agreement achieved 
by the 4p-4h description, is at the same time a blessing 
and a challenge. Surprisingly, the incorporation of all 
the other degrees of freedom must result in essentially no 
modification at all of the results. And this is not an easy 
task 0. 

The interaction SDPF.SM was originally devised for 
calculations in fixed Nhu; spaces. In a mixed calculation, 
some changes must be enforced. In the first place, we 
must retire from the pairing matrix elements of the pf- 
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number of particles in the pf-shell 

FIG. 3: Relative position of lowest states of the np-nh con- 
figurations. Diamonds represent the energies of the ground 
state and the superdeformed band-head after mixing 



shell the implicit effect of the cross shell excitations that 
now are taken explicitly into account. The amount of 
this reduction can be estimated in second order pertur- 
bation theory to be roughly the square of the off-diagonal 
matrix element (2 MeV) divided by the gap between the 
Op-Oh and the 2p-2h states (8 MeV). We model this effect 
by reducing all the pf shell T=l pairing matrix elements 
a 40%. Moreover, the closure of the ld5/2 orbit affects 
differently to the correlation energies of the np-nh con- 
figurations and we must correct the calculations for this 
limitation. To quantify the effect, we have performed cal- 
culations in the complete sd — pf space for the 2p-2h and 
4p-4h configurations (in this last case for high- J members 
of the band). The results show that the 2p-2h and 4p-4h 
configurations gain 1.5 MeV and 5.0 MeV more than the 
Op-Oh configuration when the ld5/2 orbit is opened. In 
order to simulate this effect, we have added to the pf- 
shell monopole interaction a term Se n p f + 1/2 SV n p J , 
with Se = - 0.5 MeV and SV = - 0.5 MeV. This pre- 
scription overestimates the correlation energy of the 8p- 
8h states, but it has no relevant influence in the states 
we are aiming to. Concerning the B(E2)'s, in the case of 
the 4p-4h band, the opening of the ld5/2 orbit produces 
a 12% increase of their values. 

Once the interaction has been adapted to the valence 
space, we proceed to make calculations -using the shell 
model code ANTOINE [l2l flij]- in spaces of fixed np-nh 
character. Naively one would locate these states at an 
excitation energy equal to n times the quasi-particle gap 
between the sd and the pf orbits (~7 MeV). In Fig|3 
this corresponds to a straight line with slope 7 MeV. 
This estimation misses two crucial ingredients; first, the 
presence of quadratic terms in the monopole hamiltonian 
that can lower substantially the uncorrelated energy of 
these configurations, and second, the gain in energy of the 
maximally correlated states within these configurations. 
To measure the size of the first contribution, we have 
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FIG. 4: The superdeformed band in 36 Ar; E 7 's, comparison 
with the results of the mixed calculation 

computed the energies of the lowest Slater Determinant 
for each np-nh space. We have plotted them in Fig. [3] 
Notice that the behavior is not linear; in particular the 
4p-4h configuration does not lay 28 MeV above the Op- 
Oh, but just 12 MeV (see 0] for a fully worked example 
of this mechanism) . The next step is to diagonalize the 
full interaction in the np-nh spaces to incorporate the 
correlations at fixed configuration. The results are the 
squares in Fig. [3] and the correlation energy can be de- 
fined as the difference between circles and squares in the 
figure. The key point is that these correlation energies are 
very large. In the 4p-4h space, they are large enough to 
produce a super-deformed rotor whose correlation energy 
exceeds that of the spherical ground state by 15 MeV, a 
huge effect that represents 5% of the total binding en- 
ergy. This suffices to bring the 4p-4h band-head close 
to degeneracy with the Op-Oh state. It is important to 
notice that the combined effect of the monopole hamilto- 
nian and the correlation energy favors clearly the 4p-4h 
super-deformed band which becomes the lowest config- 
uration above the Op-Oh ground state. Both the 2p-2h 
and the 6p-6h band-heads are higher in energy. This has 
important consequences for the mechanisms of mixing. 

The final step consist in the diagonalization of the in- 
teraction in the full valence space. As a consequence of 
the mixing with the 2p-2h states, the Op-Oh dominated 
(70%) ground state gains 5 MeV, while the 4p-4h domi- 
nated (70%) band-head of the super-deformed band gains 
2.5 MeV (see Fig. 0J, to finish close to its experimental 
location (4329 keV exp vs. 4319 keV th.). The relative 
placement of the configurations at fixed np-nh and the 
intrinsic structure of the np-nh states produce the fol- 
lowing very interesting effect: The 4p-4h super-deformed 
states mix mostly with the 6p-6h ones. A larger mixing 
with the 2p-2h states would have prompted an earlier and 
more intense decay-out of the super-deformed band. The 
Op-Oh and 2p-2h states are of single particle (or spherical 
nature) and the cross shell pairing interaction -i.e. the 




J 

FIG. 5: The superdeformed band in 36 Ar; B(E2)'s, exp. vs. 
mixed results. 

scattering of a nuclear Cooper pair from the sd to the pf 
shell- mixes them very efficiently. The lower 4p-4h and 
6p-6h states are both super-deformed and pairing mixes 
them too. What turns out to be severely hindered is 
the mixing of the spherical 2p-2h states with the super- 
deformed 4p-4h states, as semi-classical arguments sug- 
gest. In addition, the states of the super-deformed band 
cannot decay to those of the 6p-6h band because they lay 
at higher energies. The decay out of the super-deformed 
band must then proceed via its small Op-Oh and 2p-2h 
components. 

In Fig. 0] we compare the 7 energies along the SD band 
with the experiment. Even if it seemed difficult to im- 
prove upon the agreement of the 4p-4h calculation, the 
mixed calculation succeeds. Now the backbending is bet- 
ter reproduced while the moment of inertia is still cor- 
rect. For the first, the interaction of the SD 10 + with 
an isolated 2p-2h 10 + does the job. For the second, the 
off-diagonal pairing, active in the mixed calculation, com- 
pensates the pairing reduction of the effective interaction 
(as it should). In the range of sensible values of the pair- 
ing interaction there is surely room to absorb the hypo- 
thetical modifications of the moment of inertia of the SD 
band due to the opening of the ld5/2 orbit. 

In Fig. we compare the experimental and theoreti- 
cal B(E2)'s through the SD band. The accord is excel- 
lent and would become perfect should we increase the 
calculated values by 12%, the quantity that we have ob- 
tained as the effect of the opening of the ld5/2 in our 
fragmentary 4p-4h calculations. The reconstruction of 
the 4p-4h results is due to the fact that the 4p-4h SD 
states mix with 6p-6h states that are in turn almost as 
deformed as them. A larger mixing with the (spherical) 
2p-2h states would have produced an unwanted reduction 
of the B(E2)'s. 

In table [I] we have gathered the experimental results 
obtained in 0] for the transitions connecting states of 
the SD band with other states of Op-Oh or 2p-2h char- 
acter. The decay-out of the band occurs through the 
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TABLE I: Out-band transitions in the SD band of 36 Ar 
(B(E2)'s in e 2 fm 4 and energies in keV) 



E 7 B(E2) 
experiment theory experiment theory 



2+ — 




4950 


4846 


4.6(23) 


4.0 






4166 


3917 


2.5(4) 


1.2 






1697 


946 


19.2(30) 


18.4 




4f 


3552 


2787 


5.3(8) 


0.25 


10+ — 




1975 


1192 


43.6(74) 


13.1 


12+ — 


10+ 


3448 


3655 


15.0(30) 


1.5 



transitions 2g D — > 0J, Ag D — > 2{ and 6g D — > Aj, that 
share two characteristics; a very small transition strength 
and a large energy release. The calculation reproduces 
fairly well these properties, thus giving a microscopic ex- 
planation of the decay-out. At the bottom of the SD 
band, the energies of the in-band emitted 7's become 
small, and the phase space enhancement of the transi- 
tions to the low lying spherical states compensates the 
smallncss of the transition strengths. In the upper part 
of the band, the presence of a 2p-2h 10 + yrast state, pro- 
duces a weak out-band excursion that is also well given 
by the calculation. 

The mean field calculation of ref. 7] is possibly the 
closest to ours. Albeit it is limited to the low spin 
regime J<6, it incorporates the mixing of different pro- 
jected mean field solutions using the Generator Coordi- 
nate Method. As we have mentioned before, the standard 
mean field calculations do not treat properly the proton 
neutron pairing. In the cranking approximation this re- 
sults in a too large moment of inertia. Indeed, as shown 
for the 48 Cr case in ref. ^(|, the correct treatment of 
the neutron proton pairing halves the moment of inertia. 
Paradoxically.the moments of inertia of the SD bands 
calculated in |7| are too small. According to the authors 
this is a drawback of their method, related to absence of 
AK=0 admixtures, as discussed in ^t|- This spurious 
reduction of the moment of inertia overcompensates the 
effect of their bad treatment of the proton neutron pair- 
ing. The cancellation of this two -unphysical- but oppo- 
site contributions can eventually produce an unwanted 
"correct" moment of inertia. As pointed out by the au- 
thors of ref. 8], this is actually the case in their AMP- 
CGM calculation using the Gogny force. The in-band 
B(E2)'s of the SD band are clearly overestimated in Q 
while for the transitions out of the band their predictions 
are acceptable but less accurate than ours. 

To complete the description we have calculated a side 
band of negative parity with a 3~ band-head at 4178 keV, 
E 7 (5~ — 3")= 993 keV and E 7 (6" — 5~)= 2183 keV 
also reported in [JjJ. In order to describe the negative 
parity states we find that it is compulsory to open the 
ld5/2 orbit. As a full calculation is out of reach, we 
have computed the negative parity states in the space 



of the lp-lh and 3p-3h configurations, referring them to 
the ground state energy calculated in the space of the Op- 
Oh and 2p-2h configurations. The results are reasonable 
with the 3" at 5040 keV , E 7 (5" ->• 3")= 751 keV and 
E 7 (6" -> 5")= 2010 keV. Notice that 36 Ar exhibits at 
around 4 MeV a multiplet of nearly degenerate states of 
Op-Oh, lp-lh and 4p-4h nature, that the calculations are 
able to explain simultaneously. 

In conclusion, we have shown that "state of the art" 
shell model calculations in the valence space spanned by 
two major oscillator shells can provide an unified de- 
scription of many coexisting nuclear structures; spheri- 
cal states of single particle nature, negative parity states 
of octupole character and a super-deformed band built 
on 4p-4h excitations. In particular, we have shown why 
the super-deformed band appears at such a low excita- 
tion energy and how the indirect mixing of very different 
states governs the decay-out of the super-deformed band. 
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